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On 26 September 2022, NASA's Double Asteroid Redirection Test (DART) 
mission successfully impacted Dimorphos, the natural satellite of the binary 
near-Earth asteroid (65803) Didymos. Numerical simulations of the impact 
provide a means to find the surface material properties and structures of 
the target that are consistent with the observed momentum deflection 
efficiency, ejecta cone geometry and ejected mass. Our simulation that best 
matches the observations indicates that Dimorphos is weak, with a cohesive 
strength of less than a few pascals, like asteroids (162173) Ryugu and 
(101955) Bennu. We find that the bulk density of Dimorphos p, is lower than 
-2,400 kg m”? and that it has alow volume fraction of boulders (<40 vol%) 
onthe surface and in the shallow subsurface, which are consistent with data 
measured by the DART experiment. These findings suggest that Dimorphos 
isa rubble pile that might have formed through rotational mass shedding 
and reaccumulation from Didymos. Our simulations indicate that the 

DART impact caused global deformation and resurfacing of Dimorphos. 
ESA’s upcoming Hera mission may find a reshaped asteroid rather than a 


well-defined crater. 


Double Asteroid Redirection Test (DART) was a planetary defence mis- 
sionto demonstrate the feasibility of using a kinetic impactor to change 
the trajectory of an asteroid’. The impact was successful and highly 
effective, resulting ina reduction in Dimorphos’s orbital period around 
Didymos, which was initially 11 h and 55 min, by 33 + 1 min (ref. 2). The 
LICIACube Unit Key Explorer (LUKE) instrument onboard the cubesat* 
captured images of the system between 29 and 320 s after impact to 
reveal filamentary streams of ejecta and other complex patterns that 
expanded for several kilometres from the impact site*. Moreover, the 
dramatic brightening of the Didymos system by solar illumination 


of released impact ejecta was observed by ground- and space-based 
telescopes”** for many weeks after the impact. 

The obtained 33 min reduction in the binary orbital period? 
implies amomentum transfer to Dimorphos that exceeded the incident 
momentum of the DART spacecraft by a factor £ that ranges from 2.2 to 
4.9 depending on the mass of Dimorphos’, whichis not currently known 
but will be measured by the ESA Hera spacecraft in early 2027°. The B 
parameter is defined as the ratio of the target momentum increment 
after the impact to the impactor momentum, in the direction of the 
net ejectamomentum, andis related to the additional thrust from the 
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production of impact ejecta”’. B strongly depends on impact conditions 
(impact velocity and impact angle’) and target material properties, such 
as strength, porosity, bulk density and target surface structure”. 

Information about the spacecraft, the impact location and the 
impact angle are well understood for the DART impact (Extended 
Data Fig. 1and Extended Data Table 1). However, the mass and surface 
properties of Dimorphos are still ambiguous. It was not possible to 
directly measure the mass and bulk density of Dimorphos with DART 
or LICIACube. Instead, these parameters were estimated from the total 
mass of the binary system, as derived from Dimorphos’s orbit”, and 
updated volume estimates of Didymos and Dimorphos provided by 
DART (ref. 1). The estimated bulk density of Dimorphos ranges from 
Ps =1,500 to 3,300 kg m” (refs. 1,7). 

The surface material properties and subsurface structure of 
Dimorphos were also not directly measured. However, these target 
parameters are vital for understanding the impact process and for 
transforming the kinetic impactor method froma full-scale experiment 
by DART into a well-understood and reliable mitigation technique for 
planetary defence. Moreover, the material properties of Dimorphos 
provide information about the origin and evolution of the Didymos- 
Dimorphos system as well as the overall characteristics of rubble-pile 
asteroids and binary asteroid systems. 

In this work, we simulate the DART impact numerically and com- 
pare the results with observations to infer the properties of Dimorphos. 
We performed numerical simulations of the DART impact using the 
shock physics code Bern SPH (refs. 14,15) over a range of assumed 
sets of material properties and interior structures for Dimorphos. 
We represented the DART spacecraft as a low-density spherical pro- 
jectile of equivalent mass (‘Projectile’ section in Methods) while the 
impact velocity and angle were kept fixed. We simulated the asteroid’s 
response tothe DART impact for up to1 hafter the time of impact using 
the numerical approach developed in refs. 16,17 to model a late-stage, 
low-speed deformation (‘Modelling approach for the late-stage evo- 
lution’ in Methods). Bern SPH’s fast-integration scheme has been 
validated against laboratory experiments" and was recently success- 
fully applied to model the impact of the Hayabusa2 Small Carry-on 
Impactor”. Due to the relatively short timescales modelled, the rota- 
tion of Dimorphos around Didymos’s and Didymos’s gravity were not 
accounted for. 

We obtained realistic configurations of boulders for our rubble- 
pile targets from simulations of the gravitational collapse of a cloud 
of spherical particles with a predefined size-frequency distribution 
(SFD)"’. To closely replicate the topography described by ref. 1, we 
then selectively removed particles positioned near the surface. 

To explore a large possible range of boulder volume fractions 
(Oto 50 vol%), we replaced some of these boulders with matrix material 
when we built our target asteroid (Extended Data Fig. 2). We explicitly 
modelled only boulders larger than 2.5 m in diameter, and the space 
between boulders was filled with matrix material. Boulders smaller 
than 2.5 m were removed from the SFD because they were too small to 
be resolved individually. Thus, components smaller than 2.5 m were 
considered part of the matrix, which was modelled as a granular mate- 
rial with low but limited cohesion. This approach created an asteroid 
whose interior structure is like that of its surface and whose structure 
overall is consistent with a gravitational collapse. For the purpose of 
this study, other deep interior structures were not considered. 

The bulk porosity of Dimorphos results from macroporosity 
between individual rocks and boulders as well as microporosity within 
rocks. An analysis of reflectance spectra of Didymos indicates that the 
best meteorite analogues for boulders on Dimorphos are L/LL ordinary 
chondrites”? ”, which have grain densities of -3,200-3,600 kg m° and 
low microporosities of ~8-10% (ref. 23). Using the method described 
by ref. 24 (‘Macroporosity calculations for Dimorphos’ in Methods), we 
calculated macroporosities of 38 + 3% from the boulder SFD in the last 
complete image taken by DART (ref. 1) and 34 + 4% from the global SFDs 


Table 1 | Table of fixed and varied target parameters 


Fixed parameters 


Target dimensions' 177x174x116 m 
Target volume’ 0.00181km° 
Boulder SFD” See Methods 
Boulders at impact location’ See Methods 
Boulder tensile strength 10MPa 
Boulder porosity 10% 

Varied parameters 

Boulder packing 0-50vol% 


Grain density, p; 3,200 or 3,500kgm* 


Matrix porosity, do 35-65% 
Matrix cohesion, Yo 0-500Pa 
Matrix internal friction coefficient, f 0.4-0.7 


measured on Dimorphos” (Extended Data Fig. 3). The derived macr- 
oporosity was primarily determined by the boulder SFD (Extended Data 
Fig. 4) and the sphericity/roundness of the boulders (Extended Data 
Fig. 5). Importantly, the macroporosity estimation is largely independ- 
ent of the assumed minimum boulder size. In our simulations, the initial 
micropososity within boulders was fixed at 10% and the initial porosity 
of the matrix (macroporosity + microporosity) was varied between 35% 
and 65% (Extended Data Fig. 6). Both the porosity in the boulders andin 
the matrix were modelled using the P-a porosity compaction model” 
with parameters as summarized in Extended Data Table 2. 

Based on laboratory measurements of meteorite falls, the average 
tensile strength of ordinary chondrites is 24 + 11 MPa, and there is no 
statistical difference between L and LL types”. In all our simulations, 
the initial material properties of the boulders were kept the same and 
we employed the tensile strength and fracture model as described 
in refs. 14,15, with parameters corresponding to a tensile strength 
Y, = 20 MPa for a ~2 cm specimen. For the boulders on Dimorphos, 
the average Y, ~ 10 MPa (Table 1). However, based on previous impact 
studies (for example, ref. 18), we did not expect the impact outcome 
to be very sensitive to the boulder tensile strength. 

The response of the target matrix material to shear deformation 
is described by a simple pressure-dependent strength model**”®. The 
ability ofa material to resist different types of stresses is an indicator of 
its strength. Granular materials, for instance, may exhibit considerable 
shear strength due to the presence of van der Waals forces and the par- 
ticles’ inability to separate or slide over each other due to interlocking 
mechanisms” ”. Here we focus on the influence of the shear strength 
at Zero pressure, commonly knownas cohesion (or cohesive strength). 
Another important term inthe strength modelis the coefficient of inter- 
nal friction. Although this parameter cannot be directly determined, 
it is possible to relate it to the angle of repose (Methods) and bound 
the range of plausible values by making reasonable assumptions. The 
angle of repose of low-cohesion materials has been measured to be 
0 = 22° (f= 0.4) for glass beads”, 8 = 30° (f= 0.55) for quartz sand” and 
0 =35-45° (f= 0.7-0.9) for lunar regolith”. 

On Dimorphos, for a cohesion lower than Y, = 4 Pa, the impact 
occurs in the gravity-dominated regime in which crater growth is halted 
by the asteroid’s small gravity rather than its cohesion”. Therefore, 
we first model impacts into cohesionless rubble piles (with Y, = 0 Pa 
but with a coefficient of internal friction of f= 0.55 (ref. 31), which is 
equivalent to 0 = 30°). Given our other assumptions regarding mate- 
rial properties, we consider that these models of impacts into targets 
with no cohesion provide an upper limit on the possible momentum 
enhancement that can be achieved from the DART impact for a given 
asteroid mass. 
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Fig. 1| Momentum enhancement £ as derived from SPH simulations. 

a-d, Continuous black lines show the £ dependence on Dimorphos’s bulk 
density p, derived from dynamical simulations’ (lo uncertainty is shown by the 
dashed lines). Data points that cross a solid black line are consistent with the 
measured £ from the DART impact. The uncertainty for each simulated data point 
(vertical bars) was calculated from the absolute difference in £ calculated from 
two different methods (‘Momentum-enhancement calculations’ in Methods). 

a, Basa function of p; for cohesionless targets with the same dimensions as 
Dimorphos' and with boulder volume fractions ranging from 0 vol% (no boulders 
larger than 2.5 m) to 50 vol%. For fixed target volume (0.00181 km?) and fixed 
matrix porosity (Ø, = 45%), the mass and bulk density of Dimorphos vary with 

the boulder packing. b, £ as a function of p; for cohesionless targets with the 
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matrix porosity varying between 35% and 65% and two boulder packings: 21 and 
30 vol%. pg was calculated for a fixed asteroid volume and varies with matrix 
porosity and boulder packing. c, Bas a function of matrix cohesion (Yo) for the 
DART impact into targets with varying coefficient of internal friction (f= 0.40 or 
0.55), an assumed grain density p, = 3,200 kg m”° and 30 vol% boulder packing. 
The horizontal line shows £ derived from dynamical simulations for a target 
with p; = 2,160 kg m”? (minimum density consistent with the results ina and 
b).d, Basa function of Y, for the DART impact into targets with f= 0.40-0.70, 
Pz= 3,500 kg m” and 30 vol% boulder packing. The horizontal line shows £ 
derived from dynamical simulations for a target with p, = 2,360 kg m” (maximum 
density consistent with the results ina and b). 


The momentum enhancement £ was calculated using two dis- 
tinct methods. For the first method, £ was calculated by summing the 
momentum over all the SPH particles with ejection velocities larger 
than the escape velocity Ves. For a given impact, the magnitude of ejecta 
momentum in the direction of the net ejecta momentum is given by 


Pej = >> MeVe5| 0 (1) 


where m, and Vej are the mass and velocity vector of individual SPH 
particles, respectively. The p,; calculation takes the gravitational influ- 
ence of Dimorphos into account. However, it does not account for the 
gravitational influence of Didymos. The second method, described in 
ref. 34, tracks the velocity magnitude of the asteroid centre of mass 
postimpact by summing the momentum of all material that remained 
below the escape velocity after the reaccumulation of the ejecta. The 
absolute difference in £ resulting from the two calculation methods 
was used to calculate the error of our reported £ values. 


Results 

First, we varied the boulder volume fraction (the volume fraction of 
objects >2.5 min size) within the target between 0 and 50 vol%, while 
keeping the asteroid volume constant (Table 1). As a result, the mass 
and the bulk density of the asteroid varied with the boulder packing. 
Our simulations show that £ is relatively insensitive to the boulder 
volume fractions up to ~30 vol% (Fig. 1a). For boulder volume fractions 
larger than about 30 vol%, the number of boulders much larger than 


the projectile close to the impact point is high enough that the crater 
efficiency and, subsequently, £ are drastically reduced by boulder 
interlocking” and possibly armouring”, which hinder crater growth. 
These results indicate that at least the surface and shallow subsurface 
of Dimorphos havea low volume fraction of boulders larger than 2.5 m 
(less than ~40 vol%), which is consistent with the last few images sent 
by DART before impact’. 

Our simulation results for the DART impact on a cohesionless 
surface of Dimorphos also provide a means to constrain the bulk 
density of the asteroid using the measured £ (Fig. 1b), assuming a 
grain density in the range 3,200-3,500 kg m”. For a body witha fixed 
volume and grain density (Table 1), the bulk density is influenced by 
its porosity. The dominant effect of additional target porosity is a 
reduction in the bulk density of the target rather than a reduction in 
ejecta from compaction of the pore space. For example, decreasing 
the bulk density of Dimorphos will increase the overall crater size and 
decrease the asteroid’s mass (for a fixed volume) and escape velocity. 
This allows a greater total volume of ejecta to escape, but the ejecta 
mass is similar. Increasing the bulk density has the opposite effect. The 
consequence is that the total momentum of escaping ejecta measured 
constrains the target bulk density and porosity: for a cohesionless 
surface of Dimorphos, the upper bound on £ (-3.6) implies that the 
bulk density of Dimorphos is less than the current best estimate of 
the asteroid’s bulk density of 2,400 kg m” (ref. 1). Thus, Dimorphos is 
probably more porous and, therefore, may have arubble-pile structure 
throughout the whole body. 
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Fig. 2 | LICIACube images of the expanding ejecta cone compared with 
SPH simulation results. a, Image acquired by the LUKE instrument onboard 
LICIACube at a distance of -76 km, 160 s after the impact. b, Zoomed-in image 
of Dimorphos and impact ejecta. The approximate outline of the asteroid is 
shown in yellow. The ejecta curtain exhibits undulations, filamentary patterns 
and shadows. c, Bern SPH simulation of the DART impact into a cohesionless, 
rubble-pile Dimorphos-sized target (with f= 0.55 and ġo = 45%) at T=160 s. 
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Due to the limited visualization domain, only a portion of the ejecta are shown. 
d, Simulation of the impact into a rubble-pile Dimorphos-sized target with 

Yo = 500 Pa ($ = 2.26 + 0.28). e, Image acquired by LUKE at a distance of -71 km, 
178 s after the impact. f, Zoomed-in image of Dimorphos and impact ejecta. 

g, Same as c but at 7=178 s. h, Same as d but at 7=178 s. The optical depth of the 
ejecta cone is not computed for this comparison between the observations and 
simulation output. Panels a,b, e, f adapted from ref. 4, Springer Nature Limited. 


Since the surface strength of Dimorphos is poorly constrained, 
for a fixed boulder distribution (30 vol%), matrix porosity (Øo = 45%) 
and grain density (p, = 3,200 or 3,500 kg m”), we varied the matrix 
material cohesion (Y, = 0-50 Pa) and coefficient of internal friction 
(f= 0.4-0.7) (Table 1 and Methods). Several possible combinations of 
cohesion, coefficient of internal friction and bulk density could result 
in the observed deflection and account for the observed momentum 
enhancement (Fig. Ic,d). Despite this non-uniqueness, it is possible 
to bound the range of plausible values by making reasonable assump- 
tions since f= 0.4 is a lower limit for geological materials. For a target 
with f= 0.4 and @, = 45%, the cohesion on the surface of Dimorphos is 
probably lower than ~SO Pa (Fig. c,d). However, lower bulk densities 
(Ps < 2,000 kg m*) or higher matrix porosities (Ø, > 55%) would require 
higher cohesion to match the observations (Fig. 1). 


Ejecta curtain opening angle and morphology 
In our simulations of the DART impact into Dimorphos-like rubble-pile 
targets, we observed the temporal changes of the ejecta cone open- 
ing angle and studied the dependences on target properties. We 
found no notable dependences of ejecta cone opening angle on the 
friction coefficient of the targets. This finding contrasts with the 
strong dependence of the cone opening angle on the coefficient of 
internal friction that is found in simulations” of the DART impact into 
homogeneous planar targets. Our present simulations of impacts 
into Dimorphos-like rubble piles found that the ejecta cone open- 
ing angle and ejecta mass depend on target cohesion. Figures 2 and 
3 compare results for cohesionless targets (Y, = O Pa) with those for 
cohesive targets (Y, = 500 Pa). 

The ejecta plume for the cohesionless target (Fig. 2c,g) is more 
massive than that for the cohesive target (Fig. 2d,h). For both cases, 
the fastest ejecta, released shortly after the impact with velocities 


higher than a few tens of ms”, forma cone opening angle w = 90° 
(Fig. 3a,b). Such fast ejecta are influenced by the spacecraft geometry 
(for example, refs. 38,39), whichis highly simplified in this study. On 
the other hand, slower ejecta, released at late times after the impact 
(hundreds to thousands of seconds) with velocities Ves < V < 5m s7, 
form a wider ejecta cone angle of -140° (at 1 m s”, Fig. 3a) for the 
cohesionless target. For the cohesive target (Y, = 500 Pa), crater 
growth ceases about 100 s after the impact, before the crater grows 
large enough for the ejection angle to be influenced by target cur- 
vature. In this case, the maximum ejecta opening angle is ~120° 
(at 1 ms“, Fig. 3b). On the other hand, for the low-cohesion targets, 
the mass ejected at low velocities (Ves: < V < 10 m s”) greatly exceeds 
the low-velocity ejecta mass from the cohesive target (Y, = 500 Pa), 
implying a larger cratering efficiency and crater growth continuing 
to later times, resulting in a wider cone opening angle influenced by 
target curvature. 

The characteristics of the ejecta plume observed by LICIACube 
provide constraints regarding the target cohesion. At time after impact 
T=160 s, LICIACube’s LUKE acquired images show ejecta concentrated 
into rays, which cast shadows along the ejecta cone (Fig. 2a,b). At 
T=178 s, the bottom of the ejecta cone and the surface of Dimorphos 
are obscured by the shadow cast by the ejecta (Fig. 2e,f and Extended 
Data Fig. 7). The shadow observed at 7=178 s implies that crater growth 
and the release of low-speed ejecta continued to that time, consistent 
only with low-cohesion target cases, for example, ref. 40. 

Images obtained from LICIACube (ref. 4) and Hubble Space Tel- 
escope (HST) observations’ revealed a wide ejecta cone angle, esti- 
mated to be w = 115-139° (Extended Data Table 3). These observations 
were used to determine the ejecta opening angle at specific times, up 
to 3 min after the impact for LICIACube (ref. 4) and up to 8 h after the 
impact for HST°. To compare simulations results with observations of 
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Fig. 3 | Ejecta cone opening angle. a, Cone opening angle derived from 
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LICIACube data based on the opening angle at the base of the cone at T =170 s, 
resulting in ejecta velocities in the range of a few cm s” to a few tens of m s7 
(ref. 4). œ =125 + 10° was measured from HST data for ejecta in the range 
1-10 m s” (ref. 5). Temporal evolution measurements of the observed ejecta cone 
were not possible due to the limited observation window. 


the ejecta opening angle at a specific time after impact, we determined 
the implied ejecta velocity at the base of the visible cone in the LICIA- 
Cube images’ (Fig. 2e,f) using its distance above the surface and time 
after impact (Fig. 3). 

Overall, we found that to qualitatively reproduce the amount of 
material observed in the ejecta cone (Fig. 2c,g) as well as the observed 
cone opening angle of up to w = 139°, Dimorphos’s surface cohesion 
must not exceed ~500 Pa. From our suite of numerical simulations 
with the assumed boulder packing, matrix porosity and grain density 
summarized in Table 1, we found that the target case with f= 0.55, 
Ps = 2,200 kg m”? and Y, less than a few pascals is consistent with the 
target mechanical properties inferred from surface geology” and pro- 
duces af value (Fig. 1), excavation timescale (Fig. 2) and ejecta opening 
angle (Fig. 3) most consistent with observations. 


Deformation 

Observations from the first few hours after impact imply that more 
than 1.3-2.2 x 10’ kg of ejecta were released from the DART impact 
(equivalent to 0.3-0.5% of Dimorphos’s mass, assuming a bulk den- 
sity of 2,400 kg m°)‘. Our simulation results for weak (Y, < 50 Pa), 
Dimorphos-like targets (f= 0.55 in Fig. 1c) show that the amount of 
ejected material could be as high as 1% of Dimorphos’s mass (Fig. 4a). 
Moreover, up to 8% of Dimorphos’s mass may have been displaced 
or ejected at below the escape velocity of Dimorphos. In all impact 
scenarios simulated here, the DART impact does not produce a con- 
ventional impact crater but instead causes global deformation of the 
target (Fig. 4b). 

The outcome of the impact in terms of the postimpact target 
morphology is highly sensitive to the target cohesion. For a cohesion- 
less target, the ratio of the major to intermediate axes a/b could have 
changed fromthe reported pre-impact value of 1.02 + 0.02 (ref. 1) to as 
muchas 1.2. Such a large change in the a/b ratio is detectable with the 
highest-quality postimpact light curve data’. 

A global deformation of Dimorphos would have modified the 
gravitational field between Didymos and Dimorphos with important 
implications for its orbit. The shape change would cause an additional 
perturbation to Dimorphos’s orbit, on top of those caused by the space- 
craft momentum and ejecta recoil, and this effect can account for a 
few seconds to several minutes of the observed orbit period change 
(-33 min), depending on the magnitude of the deformation®’. Any 
deformation would alter Dimorphos’s mass distribution and affect its 
postimpact rotation state (for example, refs. 44,45). 


Implications for binary asteroid system formation 
Our numerical simulations suggest that Dimorphos is probably a 
rubble-pile asteroid with a bulk density comparable to or lower than 
that of Didymos. Our calculations based onthe observed boulder SFD 
indicate that the macroporosity estimate for the surface of Dimorphos 
(~35%) is approximately twice the value obtained for the surface of 
Ryugu, as determined through the same method”, but is compa- 
rable with the macroporosity of Itokawa***”. However, note that this 
estimate is only a rough approximation due to the limited data cur- 
rently available. The upcoming Hera mission will be able to provide 
better constrains. 

Our findings serve as crucial evidence regarding the origin of 
Dimorphos asa secondary ina double asteroid system. To maintain its 
structural stability given its rapid rotation period of 2.26 h, the primary, 
Didymos, probably requires a higher cohesive strength, estimated to be 
of the order of tens of pascals**. This level of cohesion can be attributed 
to van der Waals forces acting between the fine regolith grains” or to 
a coherent inner core®°. However, our best-fitting scenarios indicate 
that Dimorphos, the satellite of Didymos, exhibits a cohesive strength 
of less than a few pascals. This observed disparity in cohesive strength 
between Didymos and Dimorphos suggests a potential scarcity of fine 
grains within Dimorphos’s structure as well as a weak and fragmented 
internal structure. 

The material properties estimated in our study support the hypoth- 
esis that Dimorphos formed through rotationally or impact-induced 
mass shedding and subsequent reaccumulation from Didymos. The 
accretion of orbiting mass shed from Didymos occurs over a period 
of several days to years”, during which fine grains tend to escape from 
the system due to solar radiation pressure”. As a result, the accreted 
satellite, Dimorphos, has limited fines and cohesion. 

Although the mechanical properties of Dimorphos resemble those 
of Ryugu and Bennu (for example, refs. 17,53-55), these findings may 
not be applicable to single, small, S-type asteroids and specifically 
to their moons. The implications of our study may extend beyond 
Dimorphos and provide valuable insights into the formation processes 
of similar small S-type binary asteroid systems. 

Moreover, since the DART spacecraft probably caused global 
deformation of Dimorphos, we can infer that similarly formed asteroid 
moons are easily reshaped and their surfaces are relatively young“. 
Overall, the findings of this study provide valuable information for 
understanding the formation and characteristics of binary asteroids 
and will inform future exploration and asteroid deflection efforts. 
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Fig. 4 | Global deformation of Dimorphos. a, Total target mass ejected 

with speeds above v.s and total target mass displaced or ejected below Veso 
normalized by the initial target mass M; for a Dimorphos-like target with f= 0.55, 
o = 45% and cohesion levels between Y, = 0 and 50 Pa. b, Two-dimensional slices 


| 
T 
1Pa 


10 Pa 50 Pa Cohesion 


(taken at x = 0 inthe y-z plane and at z= 0 inthe x-y plane), at T ~ 1h after the 
impact. The boulder material is shown in yellow and the matrix material is shown 
in blue. The red contours show the initial target profile before the impact. 

The black arrows show the impact direction. 


Methods 

Macroporosity calculations for Dimorphos 

The bulk porosity of Dimorphos results from macroporosity between 
individual rocks and microporosity within rocks. The macroporosity 
ofagranular assembly is affected by the particle SFD, interparticle fric- 
tion and the packing history®’. Numerical simulations of gravitational 
aggregation?”*°*°, whichis one of the probable formation mechanisms 
for Dimorphos-like secondaries”, can be used to estimate the initial 
packing of the boulders within the asteroid. 

We use the method described by ref. 24 to calculate the surface 
macroporosity from the SFD of boulders at the impact location’ and the 
global boulder SFD onthe illuminated side of Dimorphos” (Extended 
Data Fig. 4a). We used a cumulative Weibull (Rosin-Rammler) distribu- 
tion” to represent the two SFDs: 

N(D) = Ny exp(—3(D/A)*)/B, (2) 
where D is the boulder size and M, A and £ are constants (Extended 
Data Fig. 4a). We found that the fitting constants N,, A and £ vary with 
the assumed minimum boulder diameter D min as seen in Extended 
Data Fig. 4b,c,d. 

We used Monte Carlo code to generate a Gaussian distribution 
for the particle bulk density (including the approximate micropo- 
rosity), the initial boulder packing, the boulder sphericity and the 
boulder roundness. For the bulk density calculations, we assumed a 
grain density of 3,500 +100 kg m”, with microporosity -10% (ref. 23), 
giving a mean value of 3,220 kg m” (Extended Data Fig. 5a). Based on 
the pkdgrav simulations described in‘Numerical model’, we considered 
an initial packing with a mean value of 39% (Extended Data Fig. 5b). 
Preliminary estimates of the sphericity and the roundness of boulders 
were derived from the last complete image taken by the DRACO camera 
onboard DART”, giving mean values of 0.74 and 0.56, respectively 
(Extended data Fig. 5c,d). Other observational output constraints are 
given in Extended Data Table 3. 


Numerical model 

In this work, we used the Bern SPH shock physics code". The 
code was originally developed by refs. 60,61 to model the collisional 
fragmentation of rocky bodies. It was later parallelized® and further 
extended by refs. 14,15 to model porous and granular materials. The 
most recent version of the code includes a tensile fracture model”, a 
porosity model based on the P-a model, pressure-dependent strength 
models” and self-gravity. The Bern SPH code has been validated ina 
number of studies (for example, refs. 63,15) and benchmarks against 
other codes (for example, ref. 64). 


We modelled the DART impact and DART-like impacts into 
rubble-pile ellipsoidal targets composed of different distributions of 
spherical boulders embedded into a matrix material (Extended Data 
Fig. 7). We tracked the evolution of the ejecta and of the target for up 
tolhafter theimpact. Due to the relatively short timescales modelled, 
the rotation of Dimorphos around Didymos and Didymos’s gravity were 
not accounted for. To obtain realistic configurations of boulders, we 
used the rubble-pile model with Dimorphos’s surface boulder SFD, 
which was generated from the gravitational accumulation simulations 
using the N-body code pkdgrav (ref. 65). 


Projectile. The DART spacecraft bus was approximately 1.2 x 1.3 x 1.3m 
in size, with structures extending to approximately 1.8 x 1.9 x 2.6m 
(ref. 1). Additionally, the spacecraft featured two large solar 
arrays, each measuring 8.5 m in length. The spacecraft weighed 
579.4 + 0.7 kg at impact’. Studies of the influence of the spacecraft 
geometry oncratering” show that a simple solid representation under- 
estimates the surface area of the impact compared to a -20 m complex 
projectile. Moreover, the penetration depth of a solid sphere is larger 
than for the spacecraft itself. However, ref. 38 showed that the projectile 
geometry affects only the very early, fast ejecta generated from within 
the coupling zone. Therefore, due to resolution constraints, here we 
simplify the projectile geometry. To match the penetration depth 
of the spacecraft bus, we modelled the projectile as an under-dense 
aluminium sphere, as opposed to a solid sphere (radius r= 0.52 m 
and bulk density 6 = 1,000 kg m”) with the equivalent mass of the 
DART spacecraft. For the large cratering efficiencies studied here, 
for which the crater cavity grows many times larger than the DART 
spacecraft, we considered that a spherical projectile is a reasonable 
approximation for a DART-like impact. 


Material model. We modelled both the boulders and the matrix mate- 
rial using the Tillotson equation of state for basalt’, with modified 
initial grain densities of p, = 3,200 and 3,500 kg m”. 

For the response of the matrix material to shear deformation, 
we applied a simple pressure-dependent strength model, typical of 
geological materials*”°, with asymptotes to a certain shear strength 
at high pressures. The Lundborg strength model describes the yield 
strengthas: 


Y, n fP 
° T+fP/Vam — Yo)’ 


(3) 


where P is pressure, f is the coefficient of internal friction and 
Yam is the limiting strength at high pressure. For the weak asteroid 
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materials considered in this study, we used a constant cohesion Y, with 
astrain-based weakening model that prevents artificial clumping (like 
the approach used in ref. 67). Our model uses a linear relation between 
cohesion Y, and total strain €,,,, and it is assumed that for €,,, > 1, cohe- 
sion is lost. When modelling the matrix material, the tensile strength 
was defined by extrapolating the yield strength (versus pressure) curve 
(equation (3)) to intersect the pressure axis. Additionally, we limited the 
maximum negative pressure to Prin 2 Yo. To compute the coefficient of 
internal friction ffrom the measured angle of response @ (refs. 32,68), 
we followed the relation described in ref. 64: 


_ 2y2sin() 


We varied the cohesion Y, between O and 50 Pa and the coefficient of 
internal friction fbetween 0.4 and 0.7. 


(4) 


Tensile strength of boulders. Based on laboratory measurements of 
meteorite falls, the average tensile strength of ordinary chondrites” is 
24 +11 MPa, with no statistical difference between the L and LL types”. 
However, measurements of individual specimens span from a few 
megapascals up to -100 MPa. In all our simulations, the initial material 
properties of the boulders were kept the same and we employed the 
tensile strength and fracture model as described in refs. 14,15, with 
parameters corresponding to a tensile strength Y, ~ 20 MPafora~2cm 
specimen. For the boulders on Dimorphos, the average Y+ ~ 10 MPa. 


Porosity model. The initial porosity of the boulders was fixed at 
10% (refs. 70,71), whereas the porosity of the matrix (which includes 
boulders smaller than 2.5 m) was varied between 35% (25% macr- 
oporosity + 10% microporosity) and 65% (55% macroporosity + 10% 
microporosity). Both were modelled using the P versus «& porosity 
compaction model””. The full description of the P versus a model 
implemented in the Bern SPH code is given by ref. 14. Here we used a 
simplified version of the Pversus a model witha single power-law slope 
defined by the solid pressure P, elastic pressure P, exponent n, initial 
distension a and distension at the transition from the elastic regime a,: 


1, if P, < P. 


a(P) = (5) 


p,-P \" : 
(Ae — v(; = ) +1, otherwise. 


We assumed that a, = a. The input parameters from the matrix 
and for the boulders are summarized in Extended Data Table 2. Without 
crush-curve measurements of the surface material, our input param- 
eters were informed by experimental quasi-static crush curves of 
Earth analogues and lunar regolith. We defined crush curves that have 
crushing strengths higher than the crushing strengths for sand” but 
lower than the crushing strength for lunar regolith’’””. Extended Data 
Fig. 6 compares our P versus a crush curves with the experimental 
quasi-static crush curves for lunar dust”, gypsum (50% porosity)”, 
Lane Mountain No. 70 sand (44% porosity)” and Ottawa sand (30-45% 
porosity)”. 


Modelling approach for the late-stage evolution. Due to the vastly 
different timescales of the shock-wave propagation and crater forma- 
tion, itis difficult to numerically model the entire crater formation on 
small, weak asteroids (~-100-1,000 m). To ensure numerical stability, 
the maximum time step in a shock physics code is limited by the Cou- 
rant criterion. In SPH, the Courant criterion requires that the time 
step dtis smaller than the simulation resolution divided by the speed 
of sound in the target, c,. To model DART-like impacts on Dimorphos, 
the maximum time step must therefore be At < resolution/c, = 10s. 
On the other hand, the crater formation and ejecta reaccumulation 
time in the gravity regime can last up to a few hours. Here we applied a 
transition to a low-speed medium inthe shock physics code calculation. 


At atime transition after the initial shock has passed, movement of 
the target material is governed only by low velocities, Umateriat < Cs, and 
we can switch to a low-speed medium (using a fast time integration 
scheme). In this step, we applied a simplified Tillotson equation of 
state in which all energy-related terms are set to zero. The remaining 
leading term of the equation of state is governed by the bulk modulus 
P= A(p/po- 1), which also determines the magnitude of the sound 
speed. At transition = 0.5, 20 and 500 s, we used A = 359 MPa, 3.59 MPa 
and 35.9 kPa, respectively. The shear modulus was also reduced pro- 
portionally. The Bern SPH fast-integration scheme has been validated 
against laboratory experiments" and was recently successfully applied 
in modelling the impact of the Hayabusa2 Small Carry-on Impactor”. 


Rubble-pile model of Dimorphos 

To obtain realistic configurations of boulders, we used the rubble-pile 
model with Dimorphos’s surface boulder SFD, which was generated 
from the gravitational accumulation simulations using the N-body 
code pkdgrav®. Extra particles near the surface were removed to more 
faithfully represent the target topography reported by ref. 1. To explore 
alarge possible range of boulder mass fractions, we removed some of 
these boulders from the pkdgrav output when we built our SPH models. 
We defined four different boulder distributions, with ~21, 30, 44 and 
50% of the target volume occupied by boulders (Extended Data Fig. 2c). 
Boulders smaller than Rmin = 1.25 m (2.5 m in diameter) were removed 
from the SFD due to being too small to be resolved individually. We 
assumed that they were part of the matrix material that was used to 
fill the voids between the larger boulders. 


Impact site and impact angle. DART impacted Dimorphos at 
8.84 + 0.45° S, 264.30 + 0.47° E, and the impact angle was 17 + 7° from 
the surface normal’. The impact was on a boulder-covered terrain, 
which was mapped using the global digital terrain model from ref. 1 
(Extended Data Fig. 1). Using Delaunay triangulation, convex hulls 
were created for the three largest boulders at the site: a 6.5 m boulder 
(Atabaque Saxum), a 6.1 m boulder (Bodhran Saxum) and a smaller, 4 m 
boulder (Caccavella Saxum). These boulders were incorporated into 
all SPH simulations. The solar arrays made contact with Dimorphos’s 
surface before the spacecraft bus; however, the majority ofthe impact 
energy was transferred through the bus, which made up most of the 
spacecraft’s mass at the time of impact’. 


Momentum-enhancement calculations 
Ina high-velocity impact event, the change in momentum of the aster- 
oid AP is amplified by the momentum of impact ejecta that escapes the 
gravitational attraction of the target body". The DART impact changed 
the orbital period around Dimorphos by -33 min and was directly meas- 
ured from ground-based observations’. This change in orbital period 
corresponds to a change in the orbital velocity of about 2.7 mm s”. 
The momentum transferred to the target body relative to 
the incident momentum of the spacecraft, quantified by the 
momentum-enhancement factor £, is defined by’: 


May = mU + m(B -1)(EOE, (6) 


where M is the mass of Dimorphos, Av is the change in Dimorphos’s 
orbital velocity due to the impact, mis the mass of the DART spacecraft, 
U is DART’s velocity relative to Dimorphos’s at impact and £ is the net 
ejectamomentum direction. 

The full equation needed to determine the momentum- 
enhancement factor f from observations is derived in ref. 9 and further 
discussed in ref. 7. We define £ as the total momentum change of the 
asteroid APdivided by the magnitude of the impactor momentum mU: 


B = AP/(mU) =1+ an (7) 
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where p,; is the magnitude of the vector sum of the momentum of 
the ejecta. 

The value of B depends on both the target material properties and 
theimpact conditions (for example, refs. 11,16,34,38,79,80). A value of 
£=1implies that the impact ejecta made a negligible contribution to 
the total momentum change, whereas £ > 2 means that the momentum 
contribution from the impact ejecta was larger than the momentum 
imparted by the impactor directly. 

To affirm the validity of our simulation results, we calculated £ 
using two distinct methods. For the first method, £ was calculated 
by summing over all the SPH particles with ejection velocities larger 
than the escape velocity Vese- Vese Was Calculated by summing over all 
the SPH particles in the simulation. However, we did not consider the 
gravitational influence of Didymos. 

For a given impact, the ejecta momentum is given by 


, (8) 


Pej = 2 MeVej 


where m, and v,; are the mass and velocity vector of individual SPH 
particles, respectively. Note that in some of the tests and analysis 
presented inthe Methods, Bis computed with respect to the along-track 
direction (y) for simplicity. In the impact simulations presented here, 
we were able to track the impact ejecta, and we performed the £ calcula- 
tion after long periods (up to 1h after the impact). Therefore, 
the p; calculation takes the gravitational influence of Dimorphos 
into account. 

Thesecond method, described in ref. 34, tracks the velocity of the 
asteroid centre of mass postimpact by summing the momentum ofall 
material that remained below the escape velocity after the reaccumula- 
tion of the ejecta. The absolute difference in £ resulting from the two 
calculation methods was used to calculate the error of our reported $ 
value for each simulation. 


Resolution tests 

So that we could model numerically the very long timescales 
required to see the impact effects, most of the SPH simulations pre- 
sented here had a limited spatial resolution of 5.6 x 10° SPH particles. 
Due to the low resolution employed here, the very fast (v/U > 10") 
ejecta was under-resolved, which may cause an overestimation or 
underestimation of the total ejected mass and ejecta momentum at 
high speeds. 

We performed resolution tests with impacts into a cohesion- 
less target with 30 vol% boulders, f= 0.55 and ġo = 45% using two 
resolutions: 5.6 x 10° SPH particles and 14.5 x 10° SPH particles. The 
low-resolution runs overestimated the amount of mass ejected at 
high velocities, M(v/U> 10°), by a few percent (Supplementary Fig. 
la). In addition, the cumulative ejecta momentum at speeds lower 
than v/U < 1e3 was overestimated by about 6% in the low-resolution 
run (Supplementary Fig. 1b). 

To further study the effects of resolution, we also computed 
the data from impacts into cohesionless targets with 30 vol% boul- 
ders, f= 0.55 and @, = 45% but where only ~15% of the target (in the 
y direction) was modelled. These targets with a reduced domain 
were modelled with 5 x 10° SPH particles and 10 x 10° SPH particles 
(Supplementary Fig. 2). We see no clear difference compared to 
Supplementary Fig. 1. 


Parameter studies 

Target curvature and impact location. For a fixed homogeneous, 
ellipsoidal target (Extended data Fig. 2, with p; = 3,500 kg m”, Y¥,=0, 
f=0.55 and @, = 45%), we considered DART-like impacts at different 
locations. Location Lis atx =0 mandy = 0 min the middle of thea and 
b axes of the asteroid. The impact was vertical. Location 2 is atx=0 m 
and z=0 min the middle of the a and c axes. The impact was vertical. 
Location 3 is the DART impact location at x = -8.31 m and z= -12.69 m 


(ref. 1). Location 4 is the DART impact location with the three boulders 
at the impact site (Extended Data Fig. 1). Both locations 3 and 4 used 
the DART spacecraft velocity vector as described in ref. 1 (the impact 
was modestly oblique). 

Supplementary Fig. 3a shows the mass-velocity distribution of 
the ejecta from impacts at locations 1-4. Supplementary Fig. 3b shows 
the cumulative ejected momentum as a function of normalized vertical 
ejection velocity (in the +z direction for location 1andin the -y direction 
for locations 2-4). For this target scenario, a vertical DART-like impact 
at location 1 produced £, = 3.82 and the impact at location 2 produced 
B,~ 3.74. At location 3, £, = 3.14, which means that the impact location 
and the impact angle reduced pajp by about 18%. At location 4, which 
has the same geometry as location 3 but has the three boulders at the 
impact site, £, = 3.06. These results suggest that the three boulders at 
the impact location had a minimal effect on £, at least for our assumed 
impactor shape. 


Boulder packing. For a fixed target shape (Extended Data Fig. 2) and 
material properties (p, = 3,500 kg m”, Y, = 0 Pa, f= 0.55 and Øo = 45%), 
we investigated the effects of the volume fraction of boulders larger 
than 2.5 m within the target on the ejecta mass-velocity distribution 
(Supplementary Fig. 4a) and the momentum-velocity distribution 
(Supplementary Fig. 4b). The boulders were modelled as spheres. 
We simulated the DART impact into targets with five different boul- 
der volume fractions: 0, 21, 30, 44 and 50 vol%. All runs included the 
three boulders at the impact site, including the target with O vol% 
boulders. We found that the DART impact into targets with a boulder 
packing of 21 and 30 vol% have almost the same ejected mass and 
ejecta momentum compared to a homogeneous, matrix-only target 
scenario (within 10%). Moreover, impacts into targets with a boulder 
packing higher than -44% show a clear reduction in the ejecta mass and 
ejecta momentum. For a target with boulder packing of 50 vol%, the 
cumulative ejecta momentum was reduced by about 40% compared 
toasimilar impact into ahomogeneous, matrix-only target. We attrib- 
ute this reduction in ejecta to a lower cratering efficiency caused by 
boulder armouring” and interlocking”. However, the DART impact 
excavated material only from the upper approximately 15 m, and our 
models cannot rule out the possibility that the interior of Dimorphos 
may be more densely packed. 

For the same impact scenario, we varied the arrangement of the 
boulders while keeping the volume boulder packing constant at 44 
and 50 vol%. The SFD distribution and the boulders at the impact site 
were the same. Our results show that the random nature of the boulder 
distribution causes a variation in the cumulative ejecta momentum of 
up to 10% (Supplementary Fig. 5). 


Grain density. Supplementary Fig. 6 shows the momentum-velocity 
distributions from the DART impact into Dimorphos-like targets with 
two grain densities, p, = 3.2and3.5 kg m”, for targets with no boulders 
(0 vol%, Supplementary Fig. 6a) and targets with 30 vol% boulder pack- 
ing (Supplementary Fig. 6b). For the DART impact scenario, the grain 
density influences only the low-velocity ejecta (v/U < 10). 


Cohesion and coefficient of internal friction. The normalized ejected 
mass at speeds greater than vasa function of the normalized ejection 
speed v/Uis shown for the DART impact into targets with varying cohe- 
sion (between 0 and 50 Pa; Supplementary Fig. 7a) and for targets witha 
varying coefficient of internal friction (between 0.40 and 0.70; Supple- 
mentary Fig. 8a). Supplementary Figs. 7b and 8b show the cumulative 
ejected momentum as a function of ejection speed. 

For the range of target cohesions investigated here, only the ejecta 
slower than -1 ms “was affected, which means that the effects are seen 
in terms of Band the amount of target deformation but notin terms of 
the morphology of ejecta at the time of the LICIACube observations 
(Fig. 2). 
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Ejecta mass and cone orientation 
Supplementary Fig. 9 shows the mass of ejecta that escapes Dimorphos 
(U > Vesce) and contributes towards the £ calculation (Fig. 1). 

Based on the analysis of the HST images taken within T + 8.5h, 
which show ejecta structures moving away from Didymos at speeds 
greater than 1ms‘, the ejecta cone centre is directed towards 
(right ascension (RA), declination (dec.)) = (141+ 8°, 25+ 6°) or 
(120 + 9°, 10 + 7°)°. The direction of the ejecta cone with projected 
velocities between a few and -500 m s” from LICIACube images’ is 
(RA, dec.) = (137 + 9°, 19 + 12°). An analysis of both HST and LICIACube 
data gives the direction of the ejecta cone as (RA, dec.) = (138°, +13°) 
with an uncertainty of 15° (ref. 7). Our analysis reveals that at T=160 s, 
the direction of the ejecta momentum vector is consistent with the 
observed cone direction within louncertainty (Fig. 3a and Supplemen- 
tary Fig. 10), which suggests that our simulation provides a reasonable 
estimate of the direction of the plume. However, our analysis of all 
ejecta, with speeds ranging from the escape velocity (V = Vese) to km s7, 
as analysed at the end of the ejecta production (T = 2,000 s), predicts 
a slightly different direction of (RA, dec.) = (150.5 + 4°, 12 + 5°) (Sup- 
plementary Fig. 10). This indicates that the direction of the observed 
segment of the ejecta plume may not necessarily reflect the true direc- 
tion of the momentum-enhancement vector, as computed from the 
entire plume’, but may be subject to a bias based on the velocity range 
at which it is observed. 


Data availability 

Additional supporting information and input data for the model 
simulations used in this work are archived on Zenodo (https://doi. 
org/10.5281/zenodo.10246671). 


Code availability 

A compiled version of the Bern SPH code as well as the necessary input 
files are available from the corresponding author upon request. The 
SPH data visualization was produced using the National Center for 
Atmospheric Research software Visualization and Analysis Platform 
for Ocean, Atmosphere, and Solar Researchers (VAPOR v.3.8.0) (https:// 
doi.org/10.5281/zenodo.7779648). 
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Extended Data Fig. 1 | Impact location. (a) Excerpt of the final full DRACO image of Dimorphos’s surface (dart 0401930049 43695 01 iof.fits), - 2s before impact’. The 
red x indicates the approximate impact location. The north pole of Dimorphos is toward the top of the figure. (b) Convex hulls of the boulders at the impact location 
extracted from the impact site DTM’. 
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Extended Data Fig. 2| Dimorphos shape and boulders SFD. (a, b) Oblate ellipsoid shape of the target, with X: 177 m; Y:174 m and Z: 116 m. The red cross at 8.84 +/-0.4 deg S, 
264.30 + / - 0.47 deg E (x ~ - 8.31, z ~ - 12.69 m) shows the impact location. (c) Boulder SFDs used in this study, with ~ 21, 30, 44 and 50% volume fractions. 
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Extended Data Fig. 3 | Macroporosity calculations from the local and the from the local and the global boulder SFD on Dimorphos, as a function of 
global boulder SFD on Dimorphos, as a function of assumed minimum assumed minimum boulder size, D min. We find 38 +/- 3% macroporosity from the 


boulder size, D nin. We find 38 +/- 3% macroporosity from the local SFD and local SFD and 34 +/- 4% from the global SFD (1-sigma uncertainty). 
34 +/- 4% from the global SFD (1-sigma uncertainty). Macroporosity calculations 
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Extended Data Fig. 4 | Boulders size-frequency distribution and Weibull parameters. (a) Cumulative number of boulders per km? as a function of size as counted in 


the final complete image from the DART spacecraft, which includes the impact site’ and on the illuminated side of Dimorphos”. (b), (c), (d) Weibull parameters derived 
by fitting Eq. (2) to the local and global SFDs. 
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Extended Data Fig. 5 | Histograms of particle bulk density, initial packing, sphericity and roundness obtained using Monte-Carlo simulations. Histograms of 


(a) Particle bulk density, (b) Initial packing, (c) Sphericity and (d) Roundness obtained using Monte-Carlo simulations. The Gaussian distributions of sphericity and 
roundness are cut at 1. 


Nature Astronomy 


Article https://doi.org/10.1038/s41550-024-02200-3 


==: Lunar dust, 33% P— a, 45% 
— LM70 sand, 44% P- a, 55% 
103 Gypsum, 50% P—a, 65% 
— P-a,10% Ott. sand, 30-45% 
= — P-a, 35% 
© 
a 
re 
a 10 
= 
5 
wn 
Ww 
foD] 
pas 
A 
10? 
10° 


1.00 1.25 1.50 1.75 2.00 2.25 2.50 2.75 3.00 


Distension, @ 
Extended Data Fig. 6 | Crush curves. Crush curves (distension, alpha, versus pressure) used in Bern SPH for boulders (10% porosity) and for the matrix material 
(35-65% porosity), compared with experimentally derived quasi-static crush curves for lunar dust”, gypsum”, Lane Mtn. #70 sand” and Ottawa sand (shaded area)”. 
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Extended Data Fig. 7 | Stereoscopic pair from LICIACube images of the 
Didymos system and from Bern SPH simulations at T ~ 178s. These 
stereoscopic pairs can be viewed in “parallel” stereo mode simply by relaxing the 
eye convergence. We suggest viewing the pair of images from a foot or so away, 
and look through the screen to infinity, allowing the two images to float across 
each other. Where the two central pictures exactly overlap, the “fused” 3-D image 


is to be found; all that is then necessary is to gently adjust the focus of the eyes, 
while the convergence remains relaxed, to obtain a clear stereoscopic image. 
This technique is called ‘Free Viewing’ of stereo pairs. For amore authentic stereo 
effect, use a Brewster format stereoscope - The London Stereoscopic Company 
OWL or similar. Detailed instructions may be found at LondonStereo.com. 
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Extended Data Table 1| DART input constraints 


Source Parameter Values Reference 


Impact conditions 


Spacecraft mass 579.4+0.7 kg [1] 
Impact speed 6.1449 + 0.0003 km/s [1] 
Impact velocity (-1.06, 5.96, 1.03) km/s [1] 
Impact angle 17 +7 from the surface normal [1] 
Impact site lat /lon 8.84 + 0.45° S, 264.30+0.47° E [1] 
Dimorphos 
characterisation 
Extent of Dimorphos 177 +2; 174 +4; 1162 m [1] 
Volume 0.00181 + 10% km? [1] 
Angle of friction from surface slope <30° [41] 
Global boulder SFD Extended data Fig. 4 [19] 
Didymos/Dimorphos spectral type S-type [22, 76] 
Didymos/Dimorphos meteorite analogues L/LL chondrites [20, 21] 
L/LL chondrites grain density 3200-3600 kg/m? [23, 77] 
L/LL chondrites microporosity 8-10% [23, 78] 
Boulder strength 10 MPa see Methods 
Global macroporosity 34 + 4% see Methods 
Impact site 
characterisation 
Local boulders SFD Extended data Fig. 4 [1] 
Boulders sphericity 0.74 +0.25 [79] 
Boulders roundness 0.56 + 0.44 [79] 
Macroporosity 38 +3% see Methods 


Observations and DART input constraints. 
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Extended Data Table 2 | Material model parameters for impact simulations into Dimorphos analogues 


Description 


Material 


Impactor 


Aluminium 


Tillotson? 


Matrix material Boulders 


Basalt 


Tillotson? 


Basalt 


Tillotson? 


Equation of state 


Grain density, pg (kg/ m?) 
Strength model 


LUND strength parameters® 
Damage strength at zero pressure, Yo (Pa) 
Strength at infinite pressure, Yjn¢ (GPa) 


1000 
von Mises 


0.34 


3200/3500 
LUND* 


0-500 
3.5 


3200/3500 
LUND 


1x 10° 
3.0 
0.8 


Internal friction coefficient (damaged), f 


Porosity model parameters (P — a)? 
Initial porosity, do 
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Extended Data Table 3 | Observational constraints 


Source Parameter Values Reference 

Ejecta plume 

morphology 
Cone opening angle from HST, w 125° + 10° [5] 
Cone opening angle LICIA, w 135° + 4° [4] 
Cone direction HST solution 1 (RA, DEC) 1414+8°, 25+6° 5] 
Cone direction HST solution 2 (RA, DEC) 120+9°, 10 7° [5] 
Cone direction LICIACube (RA, DEC) 137 +9°, 19 12° [4] 
Cone direction HST & LICIACube (RA, DEC) 138 15°, 13 15° [7] 
Ejecta mass, Mej more than 1.3-2.2x10" kg [6] 

Deflection 

efficiency 
Pre-impact orbit period (Ppre) 11 hours 55 minutes [2] 
Post-impact orbit period (Ppost) 11 hours 23 minutes [2] 
Period change, AP 33.0 + 1.0 (30) minutes 2] 
Along-track velocity change, Avr -2.7 +0.10 (lo) mm/s [7] 

] 


Momentum enhancement, 8 (3.61 + 0.20) x LOUTI [7 


Observational constrains. Here we provide a summary of the constraints derived from observations, the DRACO camera and meteorite analogues. (HST = Hubble Space Telescope 
observations’; LICIACube = LICIACube observations’). 
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